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Abstract 

The results of experimental and computer simulation studies of the yields of 
residual product nuclei in ^°^Bi thin targets irradiated by 130 MeV and 1.5 GeV pro- 
tons are presented. The yields were measured by direct high-precision 7-spectrometry. 
The 7-spectrometer resolution was 1.8 keV in the 1332 keV line. The 7-spectra were 
processed by the ASPRO code. The 7-lines were identified, and the cross sections 
defined, by the SIGMA code using the GDISP radioactive database. The process 
was monitored by the ^^Al(p,x)^^Na reaction. Results are presented for compar- 
isons between the ^'^^Bi(p,x) reaction yields obtained experimentally and simulated 
by the HETC, GNASH, LAHET, INUCL, CEM95, CASCADE, and ALICE codes. 
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1. Introduction 



Ecology has become a matter of priority in all of the developed countries. Utilization and 
management of wastes from nuclear power plants is a priority among the environmental 
protection measures. Besides, the nations that posses nuclear weapons have become 
concerned with disposal or peaceful uses (conversion) of surplus weapon plutonium and 
of highly-enriched uranium. Appropriate approaches are under development in Europe, 
North America, Japan, and Russia fl], H, |^, ^ . 

In practice, two basic approaches are under study, namely, the safe long-term storage 
of long-lived radioactive wastes, which postpones the need for destroying the wastes, 
and nuclear transmutation, which turns the wastes into stable and shorter-lived nuclides. 
Some methods that combine the two approaches are also being developed. 

The concept of accelerator-driven electronuclear facilities was suggested some four 
decades ago and developed persistently ever since is considered to be a promising tech- 
nology for waste transmutation. Such facilities differ from conventional nuclear power 
facilities that are operated in critical mode supported by delayed neutrons. An acceler- 
ator provides the neutrons to sustain a chain fission reaction, so that such facilities are 
subcritical. The neutron source is a nucleon-meson cascade initiated by 0.8-2.0 GeV pro- 
tons that strike a target located inside the facility. In terms of neutron yield, nonfissible 
Pb, Bi, W, and Hg, as well as fissible U and Th, are most suitable as target materials. 
Composite targets with a suppliment of the long-lived ^^Tc and ^°Sr fission products may 
also be interesting. 

The neutron-generating target provides a neutron flux that is sufficient for effective 
transmutation to occur in a neutron-multiplying blanket that contains nuclear fuel and 
long-lived nuclear wastes. 

The accelerator-driven facilities are advantageous mainly because of their subcritical 
mode of operation, thus precluding any emergencies accompanied by an uncontrollable rise 
of power. This is primarily due to the possibility for proton accelerators to be controlled, 
or else de-energized promptly, so that the power of the facility blanket falls completely. 

Current development of neutron-generating targets for accelerator-driven nuclear fa- 
cilities is aimed at extending the data on the numerous reactions induced by high-energy 
neutrons and protons in target materials, including generation of product nuclei within 
the range of atomic numbers A from 1 up to the target atomic number plus 1. The 
generated nuclei may be both radioactive and stable. Their nuclear characteristics may 
affect the rated performances and safe operation of a facility through: 

• total activity of the target, 

• "poisoning" the target, 

• accumulation of long-lived radionuclides to be transmuted, 

• alpha-activity of the isotopes (e.g., Po) produced, 

• production of nuclides with a high vapor pressure (T, He, Hg, etc.), 

• accumulation of chemically-active nuclides, thus reducing the corrosion resistance 
of the structure materials. 
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The present-day accuracy requirements of predicting the yields of product nuclei are 
estimated to be about 30% 0. The major difficulty is that experiments with heavy target 
nuclei can only detect from 80 to 100 product nuclei, whereas the actual nuclide production 
yields are an order of magnitude higher. It has been estimated that only 300-400 of the 
long-lived nuclides are of interest, so we have to use mainly the theoretically-predicted 
reaction cross sections. The predictive power of the computational methods and of the 
respective computer simulation codes can be improved by making comparisons with the 
still scanty available experimental data 0, §. In his review of nuclide yield data, A.J. 
Koning ^ has shown that the NSR and EXFOR databases are insufficient at energies 
above 20 MeV. 

With the above approach, the accuracy requirements may be met by step by step 
improvements in the computer simulation of intranuclear processes, as well as by extending 
the experimental data base, in particular the data on the nuclide yields in thin targets 
irradiated by up to 2 GeV protons. A comprehensive comparison between the codes and 
the experimental data is expected to help determine the ranges of energies and mass 
numbers described most adequately by the codes. In the given work the energy choice 
is caused by the presence in a target proton spectrum of two basic components which 
actively participate in secondary nuclides formation: the primary one, with energy 1-2 
GeV, created directly due to protons in the accelerator beam; and a secondary one, with 
energy basically up to ~ 0.2 GeV, formed in primary proton interactions with target 
nucleus. 

2. Techniques 

Generation of residual product nuclei in proton interactions with a target nucleus is due 
to the intranuclear processes of spallation, fission, fragmentation, and evaporation of light 
nuclei and nucleons. The general form of proton interactions with a target nucleus is 



where, as usual, (p, x) designates a nuclear reaction type, T and N are chemical symbols of 
the elements, i.e. of the target nucleus and of product nuclei. At and Zt are respectively 
the mass number and the charge of a target nuclide, A and T are respectively the mass 
number and the charge of a nuclide produced in the respective nuclear reaction. 

The present work determines the independent and cumulative yields of radioactive 
nuclear reaction products. The '^'^ Al^p^x^'^Na reaction is used as monitor. 

The variations in the concentrations of two "congener" nuclides produced in a thin 
target under irradiation may be presented as 
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^ = aiNT<i> - AiiVi , 

(2) 

= a2NT<^> + z/iAiATi - X2N2 , 

with the initial conditions A'^i(O) = and A^2(0) = 0. Here, A^^i and N2 are the numbers of 
nuclei produced, ai and a2 are their yields, Ai and A2 are the decay constants, ui and z/2 
are the probabilities for nuclide 1 to decay into nuclide 2 and 3, i.e. the branching factors, 
Nt is the number of ^^^T" in an experimental sample, $ is the proton flux and t is the 
current time. In the case of a pulsed irradiation mode, then, the solution for the set of 
differential equations (|) is described by the expressions 
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where r is the duration of a single proton pulse, T is the pulse repetition period, K is the 
number of irradiation pulses, and 
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The after-irradiation decay of the nuclides produced is described by the set 
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where Ni^ and N2g are the numbers of nuclei produced as the "cooling" starts (i.e., the 
irradiation stops). 

The number of nuclei produced as the irradiation stops corresponds to their number 
as the "cooling" starts. In this case the following condition is satisfied for K irradiation 
pulses: 

N^,=N^[{K-l)T + r] , N2,=N2[{K-l)T + r] . (8) 

Instead of the numbers of nuclei in experimental samples, i.e., the nuclear concentra- 
tion, the actual experiment measures the counting rates in the total absorption peaks 
at 7-line energies Ei and E2, i.e., the intensities in the peaks that are related to the 
concentrations Ni{t) and A^2(^) as 



= N^{t)Xir]iei , S2{t) = N2{t)X2me2 , 



(9) 
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where rji and ri2 are the yields of the 7-hnes and ei and £2 are the spectrometric effective- 
ness at 7-hne energies Ei and i?2- Applying, then the condition (H) to Si{t) and 5*2 (t) ), 
we get 



S,{t) = Aoe-^'\ (10) 
S^it) = Aie~^'' + A2e-^'\ (11) 



where 



Ao = Nr^airiiEiFi , (12) 

A, = NT^air]2e2 ^ ^\ i^iF, , (13) 
A2 — Ai 

A2 = NT^ri2£2 ( (T2 - aiui , ^\ 1 ^2 • (14) 
\ A2 — Ai / 

The coefficients Aq, Ai and A2 that carry information on the cross sections ai and a2 
are determined by least squares-fitting the experimental points Su and 821 through the 
respective functions 



git) = Aoe-^^' , (15) 

fit) = Aie-^i* + A2e-^'' . (16) 

One constructs the quadratic functionals Ri and R2 

Ri = E(Si,-Aoe-^^*-)VA^?, , (17) 

i=l 

R2 = Y.{S2^-A,e-^^''-A2e-^'''y/ASl, (18) 



i=l 

where Su and 821 are the counting rates measured in the total absorption peaks of 7- 
quanta of the first and second nuclide at moment tj, ASu and AS2i are the absolute 
errors in the above counting rates and tj is the time span from the irradiation end to the 
middle of the z-th measurement interval. If Li and L2 are the numbers of experimental 
points for the first and second nuclide, respectively, then using the condition of minimizing 
the functionals, we get the following expressions to find the parameters Aq, Ai, A2 and 
their errors: 
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'I = M-^^ , (21) 
AA, = ^{M-^h , (^ = 1,2), (22) 

where 



'1 = {A,,A2}, 1=|- 



Z2. 

The matrix M and the vector of the right-hand side 2^ of the initial set of hnear equations 
are 

L2 



M^J = E (^"^''^'^■^^VA^y , (23) 

k=l 

Z. = |:(^2fce-^'*VA^2\.) , ^,J = 1,2. (24) 

k=l 

Calculating the cross sections necessitates determination of the proton flux $. With 
that purpose, an experimental sample was irradiated together with the Al monitor sample 
for which we have, by analogy with expressions (0), (^) and (|T2|) 

SNa{t) = NM^astPNamaSNae-^''-' = ^C"^^"* , (25) 

where agt is '^'^Al{p,x)'^^Na monitor reaction cross section. The parameter B is also 
determined by least squares fitting the experimental points through a dependence of the 



form (^) using formulas (|19D and (^). The number of ^'^Na nuclei produced in the 
monitor will, then, be 

NNa = = NAl<^ast^ , 

which permits the proton flux $ to be presented as 

^^N^^ (26) 



The expressions (ll^), (P^, and (|26D can be used to calculate the cumulative yield of 
the first nuclide (or its independent yield in case its precursors are absent): 

num NAlFNa , . 

At the same time, the formulas (0), (p^, ( |21| ) and (^) can be used to obtain expressions 
for calculating the cumulative yield of the first nuclide, as well as the independent and 
cumulative yields of the second nuclide: 

cum ^1 Nai A2 - Al F^a , . 

Vim^2FlNNa Nt A2 ANa 
ind _ M2 AiAi\ 1 NAlF^a 

' ~ \F, F.X^) me^NMaNrXNa ^ ' 

ind , ,. _c«m Ml , 1 NAlFjsfa 



cum ind , cum I ^ i ^ 1 ^'/ii-^iva /'Qn^ 
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Obviously, the yields calculated by formulas (^) and ( pS]) must be the same. However, 
the yield obtained by formula (^) is usually included in the final results because the 
calculation accuracy of (^) is higher. 

3. Monitor reactions 

As noted above, the ^''Al(p,x)^^Na reaction was used monitor. The 27Al(p,x)2^Na 
cross section was calculated from the approximation function 

k 

a = Y.a,E' (31) 

i=0 

recommended by V.G. Khlopin Radium Institute with the coefficients tabulated in Ta- 
ble 1 0. 

It should be noted that Gumming |[T^ was the first (in 1963) who had obtained and 
recommended to use the values of the ^^Al(p,x)^^Na reaction excitation function. At 
present, however, the use of this reaction to monitor a proton fiux is regarded as incorrect 
because the (n,a) reaction may also contribute to the yield of ^^Na, thus overestimating the 
proton fiux density $. Therefore, the ^^Al(p,x)^^Na monitor has been preferred recently. 
The situation can be seen clearly in Fig. 1 showing the excitation functions of the two 
reactions evaluated independently at the V.G. Khlopin Radium Institute and used at 
ZSR (Germany) |jTl[|. A comparison of these two plots corroborates the advantage of 
using the ^''Al(p,x)^^Na reaction as a monitor. At the same time, this monitor does 
not seem to be very attractive for us because we have to measure the short-lived (p,x) 
reaction products, which is particularly urgent in studying medium and heavy nuclei with 
reaction products form complicated decay chains. This condition restricts the irradiation 
times of the experimental samples since only long-lived nuclides are accumulated due to 
larger amont of irradiation times and hence the loading parameters of the spectometer are 
becoming worse. For the minor accumulation of ^^Na resultant from the selected short 
irradiation times to be measured within a high accuracy, we have to use low-background 
spectrometers, thus making the research much more expensive. 

One can see from Fig. 1(A) that the data used at ZSR are above the V.G. Khlopin 
Radium Institute recommended curve by ~ 10% at 100-200 MeV, by ~ 5% at 600-1600 
MeV and by ~ 12.5% at 1600-2600 MeV. Since the accuracy of the two data groups 
is 10%, they can be regarded as coincident to within errors. Therefore, the possible 
neutron background becomes very important because its occurrence during irradiation of 
experimental samples would additionally increase the systematic differences among the 
eventual values of yields obtained by different researchers with different monitor reactions. 

4. Irradiation of Experimental Samples 

The yields of residual product nuclei from 1500 MeV proton irradiation were determined 
by exposing the experimental samples to the beam extracted by slow extraction system 
from the ITEP U-10 synchrotron. The extraction is schematically presented in Fig. 2. 
The extracted beam has the form of an ellipse with ~ 25x15 mm axes, the beam intensity 
is ~ 2 ■ 10^^ proton/pulse, the pulse repetition rate ~ 16 min~^, duration of a single pulse 
~ 0.5 s. 
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The yields of residual product nuclei from 130 MeV proton irradiation were determined 
by exposing the experimental samples to the proton beam used for medical purposes. 
The extraction is schematically presented in Fig. 3. The beam was formed by a set of 
collimators. The beam at the outlet of the last collimator is of circular form of ~ 15 mm 
diameter, ~ 5 ■ 10^ proton/pulse intensity, ~ 16 min~^ pulse repetition rate, and ~ 100 
ns single pulse duration. 

Use was made of experimental samples of metallic Bi prepared by pressing metallic 
powder. Before the exposure, the experimental samples were weighed with Sartorius BP- 
61 analytical scales and were then "soldered" together with the monitor samples and 
interlayers into polyethylene envelopes. 

In each of the irradiation runs, a monitor-Al interlayer-experimental sample sandwich 
stack was placed normally to the proton beam of the respective energy. The diameters 
of experimental sample, interlayer and monitor were strictly the same (10.5 mm). The 
exposure time of a single sample was ~15-30 min. The 1500 MeV and 130 MeV proton 
fluences were ~ 2.4 ■ 10^^ and ~ 1.5 ■ 10^^, respectively. 

After the exposure, the experimental samples and the monitors were replaced into 
sealed polyethylene envelopes to prevent the gaseous reaction products from being lost. 

The irradiation envelopes were used to check the loss of the reaction products which 
may escape from the irradiated experimental samples by analyzing the composition of the 
samples. This check demonstrated the complete absence of loss. 



5. Estimation of the neutron background 

In the present work, the possible neutron background in exposing experimental samples 
to proton beams was estimated by the SSNTD techniques used to find the parameters of 



medical proton beams |]12|. With SSNTDs, we measured the secondary radiation back- 
ground composed of neutrons and protons. In the experiment, we used the ^°^Bi,^^^U and 
^^''Np targets of different fission thresholds and the glass plate detectors with collimator 
grids. The detectors and the targets were placed near the irradiated samples at different 
distances from the proton beam. In calibration, the detectors and the targets were irra- 
diated by a proton beam for a short time. Fig. 4 shows the relative distributions of the 
track numbers from the fission fragments recorded with the SSNTDs near the extracted 
1500 MeV and 130 MeV proton beams. 

The results displayed in Fig. 4 show that the secondary background is nearly isotropic 
and that the total number of fission fragments due to the secondary radiation near the 
two beams is as small as a few hundredths of one per cent of the number of the recorded 
fission fragments in the proton beams proper. This indicates that such a background 
cannot introduce any substantial distortions to the proton flux density found using the 
^'''Al(p,x)^^Na reaction excitation function. Surely, this is valid only if the data selected 
to plot the recommended dependence do not comprise the neutron background error. 



6. Measurement and processing of 7-ray spectra 

The measurement facility is a spectrometric circuit comprising a GC-2518 Ge detector, a 
1510 module (ADC, amplifier, high- voltage supply), and a S-100 base that as an integral 
part of an IBM PC emulates a multichannel analyzer. The ^°Co 1332 keV 7-line energy 
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resolution of the facility is 1.8. keV. Fig. 5 presents the measured 7-spectra of the {p,x) 
reaction products. 

The 7-spectra measured were processed by the PC IBM-realized ASPRO code . 



The spectra thus processed were combined to form an input file for the SIGMA code. 
The code plots the intensity of a selected 7-line versus time and, using its energy and 
simulated half-life, identifies the produced nuclides by the GDISP database and simulates 
their cross sections by formulas (P7|-pil|) 

The measured 7-spectra are of high intensity, especially during the early decay period. 
Therefore, any possible detection losses under high loads of the spectrometer must be 
checked out carefully. To avoid the associated systematic errors, the spectrometer was 
checked out by the method of two sources. Fig. 6 shows the variations in the peak area 
and resolution versus the spectrometer load. 

Conforming to these results, the ultimate load of the spectrometer was limited so as not 
to exceed 5% in all of the measurements. To remain within that limitation, we started 
monitoring the experimental samples at ~500 mm height and, as the load decreased, 
descended down to the ultimate height of 40 mm. The 7-spectra measured were correlated 
to each other by introducing the height factors determined by the function R{x) = a/{b + 
x)^ with constants calculated by approximating the experimental points. Fig. 7 shows a 
plot of this function. 

The independent and cumulative yields of the reaction products in formulas ( ^ - |30|) 
were calculated using the relative spectrometric effectiveness. Fig. 5 shows the ana- 
lytical dependence of the spectrometric effectiveness in the 36.5-2650.0 keV range. The 
techniques of determining the effectiveness are described in |P3|. 

In addition, we have used in our measurements also a low-background spectrometer 
with ~50 and ~100 decay/hour ultimate sensitivities of its Ge-Li detector at 1500 keV and 
200-300 keV, respectively. This spectrometer was designed for measuring the magnetic 
moment of the neutrino and has a counterpart described in [0. This spectrometer was 
used to measure 7-spectra of experimental samples 1.5 years after the samples were irradi- 
ated. Cross sections of long-lived nuclear reaction products were obtained by identifying 
and analyzing these 7-spectra. 

The corrections for 7-absorption in the substances of experimental or monitor sam- 
ples were calculated using a expression published in [0. The values of the total linear 
attenuation factor /i were taken from [|18| . 



7. Experimental Results 

Tables 2 through 5 present the results of measuring the independent and cumulative yields 
of the ground and metastable states of the ^°^Bi nuclear reaction products at proton 
energies of 1500 MeV and 130 MeV. Table 2 presents the experimental cumulative yields 
of the ground states of the reaction products. Table 3, the experimental and simulated 
independent yields of the ground states of the reaction products. Table 4, the experimental 
and simulated independent and cumulative yields of the metastable states of the reaction 
products and Table 5, the total experimental yields of the ground and metastable states 
together with the simulated yields of the ground states of the reaction products. 

The ^^Al(p,x)^^Na cross sections were calculated according to (|3T| ) and are taken to 
be 9.9±1.0 mb and 9.6±1.0 mb at 1500 MeV and 130 MeV proton energies, respectively. 
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It should be noted that the experimental values have been renormalized to the quantum 
yields in the PCNUDAT nuclear database [0] and Table of Isotopes (8th Edition) ||2^ . 
The renormalization has become necessary because the tabulated experimental yields of 
reaction products were obtained for many of the nuclides by averaging over several 7- 
lines. When varying the recommended 7-quantum yield data, therefore, their separate 
values have to be renormalized and averaged again. The complete information cannot be 
published in practice because it is too voluminous. The reaction yields obtained using the 
low-background spectrometer are labelled (*). 

Fig. 8 shows some of the most characteristic plots of the time dependences of the 
^°^Bi(p,x) reaction product decays. 

In the case of ^°^Bi irradiated by 1500 MeV protons. Fig. 8(A) shows the time de- 
pendencies of two independent ^''^Re and ^^^Pb decays at E^~290.0 keV and of the in- 
dependent ^^Tc decay together with the decay chain ^^Zr— >-^^Nb of parent nuclides at 
E7~765.8 keV. In the case of ^°^Bi irradiated by 130 MeV protons. Fig. 8(B) shows the 
time dependencies of the independent ^°^Po decay at E7=872.4 keV, of the decay chain 
^o^Bi— >^°^Pb of parent nuclides at E7=331.2 keV. From Fig. 8 it can be seen that the 
analyzing the decay curves of the nuclear reaction products may be much more compli- 
cated compared with the simplified case described above in the section 2. This is due, 
first of all, to the fact that the nuclear transition energies of the reaction products may 
prove to be so alike that the actual unresolvable 7-lines occur in the measured 7-spectra. 
Given the situations like the above, the SIGMA code realizes the feasibility of separating 
them by the least squares method and in view of the fact that the nuclides produced are 
of different half-lives. One of the SIGMA subroutines separates the nuclide decays by 
independent decay chains, while another separates an independent nuclide decay from a 
chain of parent nuclides. 



8. Measurement errors 

Analyzing the errors in the results obtained has shown that they fall mostly within = 
(10-T-35)%. The values of the errors were calculated as follows. Since the results presented 
were obtained mostly by averaging several ((TjiAaj) values calculated for different 7-lines, 
their means and their errors were calculated using formulas 



a 



where Wi = 1/Acr, 
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(33) 



The total error in the values presented in Tables 2, 3, 4, 5 was calculated making 
allowance for the monitor error as 



Aw 
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(34) 



The errors in the independent and cumulative yields of the reaction products for 
separate 7-lines, Acxj calculated by formulas (^) and (0) were determined using the 
formulas for transfer of errors. 
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As shown by the analysis, the uncertainties in the nuclear data on 7-yields and on 
monitor reaction cross sections make the major contribution to the total error. 



9. Theoretical simulation of measured products 

9.1. Technique of comparing experimental and simulated results 

The simulation results should be compared with the experimental data making allowance 
for the fact that any of the computer simulation codes can only simulate the indepen- 
dent yields of nuclear reaction products, whereas most of the experimental data are the 
cumulative yields. So, any correct comparison must involve a procedure for calculat- 
ing a simulated cumulative yield on the basis of the simulated independent yields of the 
respective precursors. 

If the production chain of n nuclides is presented as 

i O"! i 0"2 • • • i 0-n 

^ \2\ ^ ... ^ \n\ 
the cumulative yield of the last (n-th) nuclide can be calculated as 

n 

= 2^<^i U , (35) 



n 

i=l 



where k is the fraction of decays of nuclide i that go to nuclide n, which is defined using 
the branching factors Uj as 



n-l 



k = 1, i = n . 



The above representation ensures from the conceptual physical treatment of the inde- 
pendent and cumulative yields and is corroborated by the form of formulas (PS|-PU|). The 



values of the nuclear chain branching factors have been taken from . Unfortunately, 
lack of any fresher systematized data has made us use the somewhat obsolete values of 
the factors. The data are very difficult to extract from the renewed databases, and the 
extraction procedure falls outside the scope of the present research program. Surely, use 
of any newer data may have affected the accuracy of the simulated cumulative yields of 
the reaction products. At the same time, identical values are very important to use when 
comparing the experimental data and the simulation results obtained by different codes. 

To get a correct comparison among the data simulated by different codes, the sim- 
ulation results are to be renormalized to unified the cross section of the proton-nucleus 
interaction, using available experimental data values or phenomenological formulas. Here 



we used Letaw's formulae described in to obtain the unified cross section to be 1663 
mbarn at 130 MeV and 1873 mbarn at 1500 MeV. 

If a comparison event with simulation differing from experiment by not more than a 
factor of 2 is assumed here to be the coincidence criterion, then the predictive power of 
the codes can be presented via a ratio of the number of coincidences to the total number 
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of comparison events. Another parameter proposed in |p to present code accuracies is 
the mean deviation of simulation results from experimental data: 



{H) = 10 




(36) 



where ( ) designates averaging over all of the comparison events. 
9.2. Computer simulation codes 

The products of ^°^Bi(p,x)-reaction were simulated by seven different codes, namely, 
• the CEM95 Cascade-Exciton Model code l23|, 



the CASCADE cascade - evaporation - fission transport code 
the INUCL cascade - preequilibrium - evaporation - fission code p6 
the HETC cascade - evaporation transport code P7[ , 
the LAHET cascade - evaporation - fission transport code 



the GNASH code based on the Hauser-Feshbach and preequilibrium approach |2 
the ALICE code with HMS precompound approach |3 



A detailed description of the models used by these codes may be found in Refs. [p^] - 
and references given therein; therefore, only their basic assumptions will be outlined 
below. 

CEM95 uses the Monte Carlo method to simulate hadron-nucleus interactions in the 



framework of an extended version ||2J] of the Cascade- Excitom Model (CEM) of nuclear 
reactions p3[ . The CEM assumes that the reactions occur in three stages. The first stage 
is the intranuclear cascade in which primary and secondary particles can be rescattered 
several times prior to absorption by, or escape from the nucleus. The excited residual 
nucleus remaining after the emission of the cascade particles determines the particle- 
hole configuration that is the starting point for the second, preequilibrium stage of the 
reaction. The subsequent relaxation of the nuclear excitation is treated in terms of the 
exciton model of preequilibrium decay which includes the description of the equilibrium 
evaporative third stage of the reaction. 

The cascade stage of the interaction is described by the standard version of the Dubna 
intranuclear cascade model (INC) |31[. All the cascade calculations are carried out in 
a three-dimensional geometry. The nuclear matter density p(r) is described by a Fermi 
distribution with two parameters taken from the analysis of electron-nucleus scattering, 
namely 

Pi.f) = Ppir) + pnir) = po{l + exp[{r - c)/a]} , 

where c = 1.07/1^/'^ fm, A is the mass number of the target, and a = 0.545 fm. For 
simplicity, the target nucleus is divided by concentric spheres into seven zones in which 
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the nuclear density is considered to be constant. The energy spectrum of the target nu- 
cleons is estimated in the perfect Fermi gas approximation with the local Fermi energy 
T^(r) = h'^[?>T^^p{r)fl^/{2mN), where m^v is the nucleon mass. The influence of intranu- 
clear nucleons on the incoming projectile is taken into account by adding to its laboratory 
kinetic energy an effective real potential V , as well as by considering the Pauli principle 
which forbids a number of intranuclear collisions and effectively increases the mean free 
path of cascade particles inside the target. For incident nucleons V = Vj\f{r) = Tp{r) + e, 
where Tp{r) is the corresponding Fermi energy and e is the mean binding energy of the 



nucleons (e ~ 7 MeV ||3T|). For pions, in the Dubna INC one usually uses a square- 
well nuclear potential with the depth ~ 25 MeV, independently of the nucleus and 
pion energy. The interaction of the incident particle with the nucleus is approximated 
as a series of successive quasifree collisions of the fast cascade particles (A^ or vr) with 
intranuclear nucleons: 

NN^NN, NN^ttNN, NN ^ ni, ■ ■ ■ ,TriNN , 

nN ttN, TTiV ^ TTi, ■ ■ ■ , TiiN (i > 2) . (37) 

To describe these elementary collisions, one uses experimental cross sections for the free 
NN and irN interactions, simulating angular and momentum distributions of secondary 
particles using special polynomial expressions with energy-dependent coefficients ||31| and 
one takes into account the Pauli principle. Besides the elementary processes (37), the 
Dubna INC also takes into account pion absorption on nucleon pairs 

ttNN NN. (38) 

The momenta of two nucleons participating in the absorption are chosen randomly from 
the Fermi distribution, and the pion energy is distributed equally between these nucleons 
in the center-of-mass system of the pion and nucleons participating in the absorption. 
The direction of motion of the resultant nucleons in this system is taken as isotropically 
distributed in space. The effective cross section for absorption is related (but not equal) 
to the experimental cross sections for pion absorption by deuterons. 

The standard version of the Dubna INC is described in detail in the monograph 
and more briefly, in the review and in the recent book by Iljinov, Kazarnovsky, and 
Paryev [Q. A detailed comparison of the Dubna INC with the well known Bertini INC 
developed at Oak Ridge National Laboratory [^] and with the popular version developed 



at Brookhaven National Laboratory and Columbia University by Chen et al. is made 
in Ref. |3§. 

An important point of the CEM is the condition for transition from the intranuclear 
cascade stage to preequilibrium processes. In a conventional cascade-evaporation model, 
cascade nucleons are traced down to some minimal energy, the cut-off energy Tcut being 
about 7-10 MeV, below which particles are considered to be absorbed by the nucleus. In 



the CEM [^, it was proposed to relate the condition for a cascade particle to be captured 
by the nucleus to the similarity of the imaginary part of the optical potential calculated 
in the cascade model Wopt.modX''^) its experimental value Wopt.exp.i^) obtained from 
analysis of data on particle-nucleus elastic scattering. The agreement between Wopt.mod. 
and Wopt.exp. is assumed to occur when the proximity parameter 

V =1 {Wopt .mod. ^^opt.exp. .exp. I 
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becomes small enough. In CEM95, a fixed value V = 0.3 extracted from the analysis of 
experimental proton- and pion-nucleus data at low and intermediate energies is used. 

The subsequent interaction stages (preequilibrium and equilibrium) of nuclear reac- 
tions are considered in the CEM in the framework of an extension of the Modified Exci- 
ton Model (MEM) [0, At the preequilibrium stage of a reaction the CEM takes 
into account all possible nuclear transitions changing the number of excitons n with 
An = +2, —2, and 0, as well as all possible multiple subsequent emissions of n, p, d, 
t, ^He, and ^He. The corresponding system of master equations describing the behavior 
of a nucleus at the preequilibrium stage is solved by the Monte Carlo technique |^ . 



Let us note here that in the CEM the initial configuration for the preequilibrium decay 
(number of excited particles and holes, i.e. excitons no = po + ho, excitation energy Eq, 
linear momentum Po and angular momentum Lo of the nucleus) differs significantly from 
that usually postulated in exciton models. 

To be able to analyze reactions with heavy targets and to describe accurately excitation 
functions over a wider range of incident particle energy, the CEM has been extended 
recently [^. The extended version incorporates the competition between evaporation 
and fission of compound nuclei, takes into account pairing energies, considers the angular 
momentum of preequilibrium and evaporated particles and the rotational energy of excited 
nuclei, and can use more realistic nuclear level densities (with Z, N, and E* dependences 
of the level density parameter). 

The extended version of the CEM realized in the code CEM95 is described in details 



in Ref. |2J]. A detailed analysis with CEM95 of more than 600 excitation functions for 
proton induced reactions on 19 targets ranging from ^^C to ^^''Au, for incident energies 
ranging from 10 MeV to 5 GeV and a comparison to available data, to calculations using 
approximately two dozen other models, and to predictions of several phenomenological 
systematics is presented in the next paper by Mashnik et al. in this issue. A com- 
parison of many excitation functions calculated with CEM95 with predictions of several 
other codes may be found in the recent NEA/OECD document [Q. In present work, all 
calculations with CEM95 were performed using for the macroscopic fission barriers the 
Yukawa-plus-exponential modified LDM of Krappe, Nix and Sierk and Cameron's 
shell and pairing corrections for the ground state and Barashenkov et al. |]42 



cor- 



rections for the saddle-point masses for microscopic fission barriers; the third Iljinov's et 
al. systematics for the level density parameters ||4^, with shell corrections by Cameron 



et al. p3 , without taking into account the dependence of fission barriers on angular mo- 



menta, with the dependence of fission barriers on excitation energy proposed in Ref. E5 



and with a fixed value of 1.082 for the ratio aj/an for both incident energies. All values 



of other CEM parametrs are fixed and are the same as in Refs. p3| , |2i|. 



The codes included into the CASCADE complex |^ have been developed and em- 
ployed for many years in the Laboratory of Computing Technique and Automation of 
JINR, Dubna to describe interactions of hadron and nuclei with nuclei and with gaseous 
and condensed matter. 

Calculation of intranuclear cascades is performed in diffuse nucleon clouds, the space 
distribution of which is defined from experiments on electron scattering. The model takes 
into consideration the competition of evaporation and fission processes in excited nuclei 
remaining after the completion of the cascade and preequilibrium (described, according 



to 1 38 1) processes. Inelastic interaction of the nuclei appears as a superposition of two- 
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particle nucleon-nucleon and pion-nucleon collisions which can occur both in the nuclear 
overlap region and inside each nucleus separately (see the second paper in Ref. p5[). In- 
teraction of cascade particles with intranuclear nucleons occur also in the case when the 
interacting nuclei are already separated and continues till all of the cascade particles are 
either absorbed by the nuclei or escape the nuclei. Mutual interactions of the cascade 
particles are neglected. The number of nucleons in the clouds (Fermi particles) is equal to 
mass numbers of striking nuclei in the beginning of the process and gradually decreases 
due to knocking out of the nucleons by the growing shower of cascading particles (the 
so-called "trawling" effect). Therefore, the coordinates of all nucleons are simulated in the 
beginning of the collision and are redefined during the calculation procedure depending 
on the cascade development. At a starting point t, corresponding to a moment when the 
colliding nuclei come into contact, all particle collisions permitted by kinematics and the 
Pauli principle are simulated. Only one which takes place earlier than the others is chosen 
{dt = min{ti}); after that the positions of the interacting nucleons and all cascading par- 
ticles (nucleons and produced pions) are shifted to the new points corresponding to the 
moment of time t + dt. The calculation procedure is further repeated until all the cascade 
particles in the colliding nuclei are exhausted (see details in the third paper of Ref. [pS]). 
Decay of the excited nuclei is calculated basing on the evaporation model as described 
in with a fixed level density parameter which is taken to be a=A/10 [MeV~^]. In 
order to consider the internuclear cascades in material samples the codes are supplied 
with simulation of the particle transport including energy losses on ionization, definition 
of the interaction point considering the boundaries of the media with different chemical 
composition. Neutron transport is calculated basing on the 26-group system of neutron 
cross-sections. 



The cascade - preequilibrium - evaporation - fission code INUCL was developed during 



several years at ITEP, Moscow by Stepanov with co-authors |£6|. The cascade stage 
of INUCL was inspired by the standard Dubna INC but differs from it in several 
points, as: use of new, more complete than in experimental data for cross sections of 
elementary NN and ttN processes (37), the range of incident particle energies was extended 
in INUCL up to 15 GeV (in the standard Dubna INC, it is only up to about 5 GeV), use 
of a different radial distribution of nuclear pairs absorbing pions inside a nucleus (38), 
taking into account the local reduction of nuclear density ("trawling" effect) at energies 
above several GeV and several other details (see the first paper of Ref. pB|). 

The preequilibrium stage of INUCL was inspired by the MEM ^ and GEM 
but differs from the Dubna version by using another approach for the squared matrix 
elements of the transition rates and by neglecting emission of complex particles and not 
taking into account the forward peaked angular distributions of preequilibrium nucleons 
(see details in the second paper of Ref. PB[] ). 

The evaporation model in INUCL was also inspired by the Dubna version |31|, |38| , 
except the use of other level density parameters and, at low excitation energies (below 
the separation energy of nucleons or complex particles), calculation of emission of low- 
energy gammas (see details in the second paper of Ref. |l26|). 

The most interesting point of INUCL is incorporation at the evaporation stage of re- 
actions of a new, "home-grown" thermodynamical model of high energy fission (see the 
third and fourth papers in Ref. p6|). From a physical point of view, the thermodynamical 
model of fission is statistical by nature, but provides A-, Z- and energy-distributions of 
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fission fragments quite different from tlie well known Fong's statistical model of nuclear 



fission uses a own systematics for the ratio of level density parameters aj/an^ and 
seems to describe experimental data better (see details in the third and fourth papers of 
Ref. m). 



The popular transport cascade - evaporation code HETC was developed at Oak 
Ridge National Laboratory. HETC is basically an extension of the code NMTC [O 



to allow particle transport up to several hundred GeV |2^. For nonelastic collisions at 
energies below ~ 3 GeV, HETC uses Bertini INC model |Q at the cascade stages of 
reactions and Guthrie's evaporation program |Q to determine the energy and direction 
of emitted cascade nucleons and pions and evaporated n, p, d, t, ^He and a-particles, and 
the mass, charge, and recoil energy of the residual nucleus. 

As we mentioned above, a detailed comparison of the Bertini INC [0 with the Dubna 



INC [PH and with Chen's et al. INC PH| may be found in Ref. PB|. We use on our work 



the standard version of HETC [0, well known and widely used in many laboratories. 



therefore we do not consider necessary to repeat here its detailed description. 

LAHET is a Monte Carlo code for the transport and interaction of nucleons, pions 
muons, light ions, and antinucleons in complex geometry [Q; it may also be used without 



particle transport to generate particle production cross sections. LAHET is the result of 
a major effort at Los Alamos National Laboratory to develop a code system based on 
the LANL version of the HETC Monte Carlo code for the transport of nucleons, pions 
and muons, which was originally developed at Oak Ridge National Laboratory p7| , ^ . 



Due to many new features added at LANL, the code has been renamed LAHET, and the 
system of codes based on LAHET designated as the LAHET Code System (LCS) PB[ . 

LAHET can use Bertini INC |^4[ (from HETC) to describe nucleon-nucleus interac- 
tions below 3.5 GeV, and a scaling law approximation to continue the interaction energy 
to arbitrary high energies, although a reasonable upper limit is about 10 GeV. 

As an alternative to the Bertini INC, LAHET contains the INC routines from ISABEL 
code. The ISABEL INC model is an extension by Yariv and Fraenkel [BO] of the VEGAS 



code |35[ mentioned above. ISABEL has the capability of treating nucleus-nucleus inter- 
actions as well as particle-nucleus interactions. In the present work, we use both Bertini 
INC and ISABEL versions of the INC to perform calculations with LAHET. 

Let us mention here two more points between the interesting features of LAHET 
(see details in [^). First, LAHET allows one to calculate preequilibrium processes as 
an intermediate stage between the intranuclear cascade and evaporation/fission in the 
framework of a multistage multistep preequilibrium exciton model ("MPM") [^. As 
initially suggested and used in the MEM the MPM uses the Monte Carlo method 

to solve the system of master equations describing the process of equilibration of excited 
residual nucleus remaining after the cascade stage of a reaction. Nevertheless, there 
are several important differences between MPM and MEM. First, the master equations 
of the MPM is simplified as compared to the one of MEM: MPM takes into account 
only nuclear transitions with An = +2, i.e., only in the direction of equilibration, while 
the MEM considers all possible transitions An = +2, = —2, and =0, accounting for all 
possible positions of the particle-hole pairs with respect to the Fermi level (transitions 
with An = 0). Then, the master equations of the MPM does not take into account the 
angular distributions of preequilibrium particles (it should be noted, that as an option. 
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the MPM and correspondingly the code LAHET, allow to calculate angular distributions 
of preequilibrium particles by using the phenomenological parametrization of Kalbach 
(see [§T[)), while the MEM does this and the version of the MEM used in CEM95 takes 



into account the conservation of momentum and angular momentum of the nuclear system 
at both preequilibrium and evaporation stages of a reaction. Then, the Monte Carlo 
algorithms used in the MPM amd MEM for solving the corresponding master equations 
also differ. There are several other differences between the MPM and MEM as using of 
different approximations for the inverse cross sections and Coulumb barriers, for the level 
density parametrs and for the matrix elements of nuclear transitions, etc. (see details 
in ^ 0)- Note, that the interface between the intranuclear cascade model and the 



MPM and the interface between MPM and the evaporation model used in LAHET 
and the ones used in CEM95 are also different. 

Let us note also that LAHET includes as user options two models for fission induced by 
high energy iteractions: the ORNL model |B2[ , and the Rutherford Appelton Laboratory 



(RAL) model by Atchinson [Q]; the fission models are employed with the evaporation 
model. The RAL model allows fission for Z > 71, and is the default fission model in 
LAHET. The RAL model really is two models, for actinide and for subactinide fission. A 



detailed description of the models incorporated in the LCS may by found in Ref. |28| and 
is available on the World Wide Web in several documents (see p8[]). 



The GNASH code applies the Hauser-Feshbach theory to calculate the decay of 
compound nuclei in an open-ended sequence of decay chains, conserving angular momen- 
tum and parity at all stages of the reaction. Prior to equilibrium decay, the preequilibrium 
emission of fast particles is accounted for using the exciton model, for up to two preequilib- 
rium ejectiles prior to equilibration. The exciton model also includes angular momentum 
considerations to determine the spin distributions of residual nuclei after preequilibrium 
emission Transmission coefficients are calculated from the optical model, and the 

continuum (statistical) level density formalism of Ignatyuk et al. is matched onto exper- 
imental low-lying discrete levels. In this way, radionuclide production can be calculated 
to both the ground state, and isomeric states. 

The ALICE code calculates equilibrium decay with the Weisskopf-Ewing theory, and 
preequilibrium emission is determined using the new Hybrid Monte Carlo Simulation 
(HMS) model [^|. This model follows successive interactions of excited nucleons creating 



three quasi-particle excitations during the equilibration process. Importantly, the theory 
is able to account for any number of multiple preequilibrium processes in a natural way - a 
major advantage over some other preequilibrium models which only consider a maximum 
of two preequilibrium ejectiles. The model has been applied successfully at energies up 
to 400 MeV for the calculation of excitation functions 



Since most of simulation codes (except for GNASH) cannot simulate the metastable 
states of product nuclei, the respective nuclide chains used to simulate the cumulative 
yields of the product nuclei were taken to be simplified. 

The results obtained by the codes are displayed in two sets of figures, first, in Figs. 
9 and 10, which show the product mass distribution for Ep = 130 MeV and 1500 MeV, 
respectively, and, second, in Figs. 11 and 12, which present the simulation-to-experiment 
ratios versus relative mass difference between an initial target nucleus and a particular 
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product nucleus for Ep — 130 MeV and 1500 MeV, respectively. Together with the 
simulated results, Figs. 9 and 10 show experimental data obtained by summing the yields 
with respect to all the available isobar decay chains making allowance for their cumulating 
effect. 

The Ep = 1500 MeV experimental results shown in Fig. 10 are a visual demonstration 
of the different production channels, namely, due to spallation reactions for masses ranging 
from ~125 to 210 and due to fission reactions for masses ranging from ~30 to ~140. 
The range of masses from ~125 to ~140 is most probably transient, where the reaction 
products can be generated by both channels. The Ep = 130 MeV experimental results 
shown on Fig. 9 demonstrate only the production channel due to spallation reactions for 
masses ranging from ~190 to 207. The lack of experimental data for the fission product 
masses is accounted for by the poor exposure conditions of the experimental samples. 

Table 3 presents the experimental and simulated independent yields of the ground 
states of ^'^^Bi(p,x)-reaction products at Ep = 1500 MeV and 130 MeV for some of the 
nuclear reaction product yields. Table 4 presents the experimental and simulated values 
of the independent and cumulative yields of the metastable states of ^°^Bi(p,x)-reaction 
products at Ep = 1500 MeV and 130 MeV. The simulation results relate only to Ep — 
130 MeV since the GNASH code is valid for simulation of proton interactions with initial 
energy below threshold of pion productions (~ 0.2 GeV) 

Table 5 presents the total experimental values of the ^°^Bi(p,x)-reaction product 
ground and metastabe states at Ep = 1500 MeV together with the respective simulated 
ground states. 

9.3. Comparison between experiment and simulation 

Table 6 presents the information that shows the predictive power of each code for both 
energies, namely, coincidence statistics which include number of experimentally measured 
products Nexp, number of simulated products that were measured Ng, number of "co- 
incidences" between simulated and experimental values Nq and the mean deviation (H) 
of simulation results from experimental data 

Obviously, the simulation results shown on Figs. 9 and 10 do not contradict the 
experimental data if the simulated values run above the experiment and repeat their 
general trend. The explanation is that the direct precision 7-spectrometry technique used 
in this work makes it possible only to identify the radioactive products of a high yield 
that includes a significant part of the total yield of a given mass. Unlike the cumulative 
yields, the independed yields makes possible to performe the direct comparison with 
corresponding calculated values. 

As seen from Tables 3, 4 and 5 the majority of the simulated results agree well with the 
experimental values. However, in some cases the simulation results are underestimated 
significantly comparing with the experimental data, especially for short-lived products. 
There are a set of particular products whose simulated yields described by the majority 
of the used codes are below the experimental ones by several orders of magnitude. 

The following conclusions concerning the scope of applicability of the simulation codes 
can be drawn from comparing simulation and experiment results: 

• INUCL describes quite reliably spallation products with masses above A=180 and 
fission products with masses below A=120; the products in the intermediate mass 
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region A=120 to 180 are underestimated basically by a factor of 2 to 5; analysis 
of particular products shows that the code underestimates significantly, up to a 
few orders of magnitude, the short-lived products (see, e.g., the yields of ^°^Po and 

202PO), 

• HETC properly describes the spallation region, though, at Ep = 150 MeV, a too 
large underestimation for production of ^°^Pb and, especially of "^^^Tl was obtained; 
since the available version of the code does not include the fission processes, it 
cannot be used for predicting fission products, 

• CEM95 also predicts only spallation products but does it most perfectly compared 
with the other codes; a new version of the CEM realized in the code CEM97 (now 
in progress, see [^) is expected to be even more succesful in the spallation region 
and is extended to describe also production of daughter nuclides in the fission and 
fragmentation regions, 

• CASCADE describes satisfactorily the production of most nuclides in the spallation 
region but overestimates some yields of nuclides in the intermediate mass region 
from A=178 to about 195 and strongly underpredicts production of several short- 
lived nuclides (see, e.g., the yields of ^oo^jj^ igsrpj ^^le spallation region and the 
yield of ^^^Eu in the intermediate region), 

• LAHET predicts the yields of most measured nuclides within a factor of 2 both 
in the spallation and fission regions, as well as in the intermediate mass region; 
nevertheless, some too big discrepancies were obtained for "^^^Tl (at Ep = 130 MeV) 
in the spallation region, for ^^^Eu and ^"^^La in the intermediate region, and for ^^Br 
and ^^As in the fission region, 

• GNASH is the only code used here capable of calculating metastable states and has 
predicted the yields of 202m. igrmp^^ well; it also described within a factor of 
2 the production of most ground state isotopes in the spallation region, except the 
yields of ^o^Bi and ^oipb, 

• ALICE also predicted the yields of most isotopes in the spallation region within a 
factor of 2, except production of ^o^Pb, ^0071, and ^^^Tl. 

10. Conclusion 

This study is the first step in our work on non-fissible targets of interest for accelerator- 
driven facilities. Final conclusions about the predictive properties of the codes should 
thus be drawn in a following paper. At the present time, we can draw only a preliminary 
conclusion that theoretical yields predicted by different codes differ sometimes up to two 
orders of magnitude. This is a strong indication that further development of all codes is 
necessary before they can become reliable predictive tools. 
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Table 1 

Approximation coefficients for the function a — Z^f=o '^i^^ {'^ ["^^] > ^ [^^^D 



E 


100-350 MeV 


350-800 MeV 


800-2600 MeV 




k=4 


k=5 


k=4 


ao 


.16082579E+2 


.91163994E+0 


.10809485E+2 




-.98011691E+2 


.72710174E+2 


.12760849E+1 


02 


.48105521E+2 


-.22074848E+3 


-.25036816E+1 


03 


-.92886532E+2 


-25781415E+3 


.10439362E+1 


04 


.63086480E+3 


-.25781415E+3 


-.13817290E+0 


05 




.77078179E+2 
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Table 2 

Experimental cumulative yields of the ground states of ^°^Bi(p,x)-reaction products at 
Ep^ 1500 MeV and 130 MeV 







£;„=1500MeV 


£;„=130MeV 


Nucleus 


Half life 


Cross Section [mb] 


Cross Section [mbl 

L J 


*207gj 


32.2 y 


58.0 ± 7.5 




206Bi 


6.243 d 


26.2 ± 3.2 


88.1 ± 9.3 


205Bi 


15.31 d 


24.5 ± 2.9 


109 ± 23 


204Bi 


11.22 h 


25.2 ± 3.0 


97.0 ± 11.4 


203Bj 


11.73 h 


19.4 ± 2.5 


107 ± 13 


202Bi 


1.67 h 


11.0 ± 2.2 


138 ± 15 


20lBi 


108 m 




85.5 ± 11.7 


200Bj 


36 m 


8.4 ± 1.0 


61.6 ± 8.0 


203pb 


51.837 h 


44.3 ± 4.9 


230 ± 35 


201pb 


9.33 h 


29.9 ± 5.0 


139 ± 17 


200pb 


21.5 h 


28.4 ± 3.3 


122 ± 14 


199pb 


90 m 


30.7 ± 5.8 


113 ± 23 


198pb 


2.4 h 


25.0 ± 5.0 


20.0 ± 8.2 


201 


72.912 h 


38.7 ± 4.7 


166 ± 27 


200 


26.1 h 


36.2 ± 4.2 


108 ± 13 


199 


7.42 h 


32.9 ± 4.5 


101 ± 18 


198 


5.3 h 


27.4 ± 4.6 


22.0 ± 3.7 


197rp| 


2.84 h 




13.4 ± 2.4 


196rp| 


1.84 h 




18.0 ± 5.3 


194 


33.0 m 


11.4 ± 3.5 




203Hg 


46.612 d 


1.15 ± 0.12 




195Hg 


9.9 h 


12.6 ± 3.8 




192Hg 


4.85 h 


23.4 ± 3.6 




198Au 


2.696 d 


0.46 ± 0.18 




i92Au 


4.94 h 


30.6 ± 4.9 




191pt 


2.9 d 


20.2 ± 4.1 




188pt 


10.2 d 


25.1 ± 2.9 




188 J, 


41.5 h 


24.7 ± 2.8 




186 jj. 


15.8 h 


10.2 ± 1.8 




184 jj. 


3.02 h 


16.4 ± 3.7 




*1850g 


93.6 d 


29.3 ± 3.6 






22.1h 


23.4 ± 2.9 






70.0 d 


17.3 ± 2.9 






12.7 h 


26.3 ± 2.9 




i«iRe 


19.9 h 


20.5 ± 4.5 






19.7 m 


16.4 ± 2.0 




i^^Ta 


8.08 h 


22.6 ± 5.4 
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Table 2 (continued) 



Ep^l500MeV £;p=130MeV 



Nucleus 


Half life 


Cross Section [mb] Cross Section [mb] 




1.18 h 


15.6 ± 2.1 


i73Ta 


3.14 h 


16.8 ± 3.1 


i^^Ta 


36.8 m 


6.6 ± 1.1 


173Hf 


23.6 h 


25.6 ± 3.1 


*172jj£ 


1.87 y 


13.2 ± 1.5 


170Hf 


16.01 h 


18.5 ± 2.6 




1.37 y 


23.2 ± 4.4 


^^^Lu 


8.24 d 


22.4 ± 2.5 


le^Lu 


34.06 d 


16.2 ± 2.7 


166Yb 


56.7 h 


11.1 ± 1.3 




28.58 h 


9.6 ± 1.1 




8.14 h 


15.1 ± 1.8 




2.34 d 


8.2 ± 1.1 




241.6 d 


2.7 ± 0.8 




17.609 d 


7.9 ± 0.9 


149Gd 


9.4 d 


8.3 ± 0.9 




38.1 h 


6.6 ± 0.9 




48.27 d 


6.1 ± 0.7 


147Eu 


24.0 d 


9.5 ± 1.2 


i^^Eu 


4.59 d 


6.7 ± 0.7 


i^^Eu 


5.93 h 


3.8 ± 0.6 


*143pj^ 


265 d 


1.07 ± 0.35 


i39Ce 


137.66 d 


2.4 ± 0.4 


i^'^La 


40.280 h 


1.3 ± 0.2 


i^'^Ba 


12.746 h 


0.82 ± 0.11 


i27Xe 


36.4 d 


2.0 ± 0.2 


121rpg 


17.78d 


1.4 ± 0.3 




115.09 d 


0.75 ± 0.12 


ii^In 


2.83 d 


1.13 ± 0.15 


105 j^h 


35.36 h 


4.3 ± 0.5 


95Tc 


20.0 h 


2.4 ± 0.5 




2.75 d 


4.5 ± 0.6 


95 Nb 


35.02 d 


6.2 ± 0.9 


952r 


64.02 d 


1.65 ± 0.21 


892J. 


78.4 h 


2.7 ± 0.4 


88 Zr 


83.4 d 


0.22 ± 0.03 


*88y 


106.61 d 


2.7 ± 0.5 


87y 


80.3 h 


4.4 ± 0.5 


85SJ- 


64.08 d 


5.2 ± 0.8 


83Rb 


86.2 d 


4.5 ± 0.9 
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Table 3 

Experimental and simulated independent yields (in mb) of the ground states of ^°^Bi(p,x)- 
reaction products at Ep^ 1500 MeV and 130 MeV 



£;„=1500MeV 



Nucleus 


Half life 


Exp. cr. sec. 


CEM95 


LAHET 


INUCL 


HETC 


CASCADE 


ALICE 


GNASH 




8.8 d 


3.3±0.5 


1.3 


4.3 


1.2 


4.5 


2.3 






202rjiJ 


12.23 d 


4.8±0.7 


4.8 


4.6 


8.6 


7.3 


3.1 






200'pj 


26.1 h 


7.8±1.7 


6.1 


5.5 


12.1 


7.1 


4.6 






146Eu 


4.59 d 


0.59±0.12 


0.41 


5.1 


0.15 




0.07 








40.28 h 


0.45±0.17 




0.08 












*102Rh 


2.9 y 


0.69±0.12 




1.5 


2.9 










95Nb 


35.02 d 


3.8±0.5 




1.2 


5.8 










76 As 


26.32 h 


3.7±0.6 




0.80 


2.7 










74As 


17.77 d 


2.1±0.3 




1.7 


2.5 










Ep=130MeV 


207po 


5.80 h 


24.5±3.8 


34.8 


38.2 


20.6 


36.2 


45.2 


35.8 


33.6 


206po 


8.8 d 


44.6±11.9 


39.1 


47.2 


28.0 


38.2 


68.5 


58.2 


32.1 


205po 


1.66 h 


30.5±3.7 


45.9 


43.7 


25.7 


39.5 


60.2 


57.8 


21.2 


204po 


3.53 h 


43.1±5.2 


45.2 


56.3 


22.5 


47.2 


85.6 


72.8 


25.5 


203po 


36.7 m 


38.5±8.9 


47.6 


52.1 


3.2 


51.0 


66.8 


57.1 


24.6 


202po 


44.7m 


50.8±7.0 


37.5 


62.3 


0.10 


57.9 


86.3 


67.3 


36.3 


206Bi 


6.243 d 


57.7±6.5 


79.1 


71.4 


66.0 


50.4 


122.5 


113.4 


168.1 


201pb 


9.33 h 


46.0±11.9 


58.9 


28.7 


25.1 


4.7 


18.1 


10.7 


10.6 


200 rpj 


26.1 h 


15.4±8.1 


6.2 


1.9 


16.5 


0.07 


0.17 


0.46 






5.3 h 


4.0±3.0 


2.5 


1.4 


18.8 




0.47 


0.96 
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Table 4 

Experimental and simulated independent and cumulative yields of metastable states of 
209Bi(p,x)-reaction products at £;p=1500 MeV and Ep^lSO MeV 







Ep— 


:1500MeV 




Nucleus 


Half life 


Yield 


Experimental 


Calculation 








cross sec. [mbl 

L J 


(GNASH) [mb] 

V / L J 


204mpb 


67.2 m 


ind 


4.4±1.3 




202mp|^ 


3.53 h 


ind 


14.0±2.1 




197 m 


44.6 m 


cum 


8.8±2.5 






1.41 h 


ind 


15.9±3.2 




194m'ji| 


32.8 m 


ind 


4.4±0.6 




193m jjg 


11.8 h 


ind 


9.9±1.4 




183m Qg 


9.9 h 


cum 


11.8±1.8 




150ml rpjg 


3.48 h 


ind 


4.4±1.0 




lOlmj^l^ 


4.34 d 


cum 


1.8±0.2 




90mY 


3.19 h 


ind 


2.4±0.4 




82mf^b 


6.472 h 


ind 


2.5±0.4 








Ep= 


= 130MeV 




202mpjg 


3.53 h 


ind 


10.8±2.3 


11.4 


197mpj_j 


44.6 m 


cum 


13.4±2.4 


10.6 



Table 5 

Experimental total yields [mb] of the ground and metastable states of ^°^Bi(p,x)-reaction 
products at Ep=1500 MeV 



Nucleus 


Half life 


Yield 


Exp. cr. sec. 


CEM95 


LAHET 


INUCL 


HETC 


CASCADE 


203 pjg 


51.837 11 


g+ml+m2 


25.9±5.8 


13.8 


11.2 


18.8 


17.0 


10.3 


196Au 


6.183 d 


g+ml+m2 


0.80±0.13 


1.6 


0.93 


2.4 


0.22 


1.2 


i^^Au 


38.02 h 


g+ml+m2 


1.37±0.24 


2.8 


1.5 


3.6 


0.30 


1.8 


96Tc 


4.28 d 


g+m 


0.96±0.22 




1.0 


2.8 






82Br 


35.3 h 


g+m 


3.2±0.4 




0.40 


2.3 






*60Co 


5.2714 y 


g+m 


1.7±0.4 




1.2 


0.64 
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Table 6 

Statistics of simulation-to-experiment comparisons. 



E„=130MeV E„=1500MeV 



Code N|+"^=30, Ni,^p=28 N|\"^=100, N|;^p=89 





1 OTTOT* 1 n 

IclLci 111 

Fig. 11 


N^lNr, 




IcLLci 111 

Fig. 12 








A 


zo / zo 


1 7"^ 
1. / o 


A 

r\. 


4Q /fi4 








1 /98 

13/23*2 


111 
2.59 




o'±j oD 


9 Qf; 
z.yo 


CASCADE 


c 


22/28 


2.89 


c 


50/64 


1.83 


HETC 


D 


19/26 


3.22 


D 


41/62 


2.64 


LAHET(ISABEL) 
LAHET(Bertini) 


E 
F 


24/28 
23/28 


1.86 
1.85 


E 


63/89 


2.51 


ALICE 


G 


23/28 


2.43 








GNASH 


H 


8/10 


1.98 









*^ N^"*^ means number of all measured products, ^%xp means number of ground state 
measured products 

*2 with exception of ^ooRi, ^oiRi, 20231, 202pQ^ 203po 
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Figures Captions 

Fig. 1. Plots of the recommended values of the monitor reaction excitation functions (A 
- 27Al(p,x)24Na, B - 27Al(p,x)22Na). 

Fig. 2. A scheme of the extraction of the 800-2600 MeV external proton beam. 

Fig. 3. A scheme of the extraction of the 70-200 MeV external proton beam. 

Fig. 4. Relative distribution of the track numbers characterizing the background of 
secondaries near the extracted 1500 MeV and 130 MeV proton beams. 

Fig. 5. The 7-spectra of the 2'^^Bi(p,x)-rcaction products measured at proton energies 
of 1500 MeV and 130 MeV (the "cooling" times are 20 hours and 3.5 hours, respectively). 

Fig. 6. Loading characteristics of the spectrometer. 

Fig. 7. Characteristics of the spectrometer (A - effictiveness, B - height factor). 

Fig. 8. Time depenedences of the ^^^Bi (p,x)-reaction product decays at Ep=1500 MeV 
(A) and Ep=130 MeV (B). 

Fig. 9. Product mass distribution of the 2°^Bi(p,x)-reaction for Ep — 130 MeV. 

Fig. 10. Product mass distribution of the 2°^Bi(p,x)-reaction for Ep — 1500 MeV proton. 

Fig. 11. Simulation-to-experiment ratios versus relative mass difference between an 
initial target nucleus and a particular product nucleus for Ep — 130 MeV. Letters are 
explained in Table 6. 

Fig. 12. Simulation-to-experiment ratios versus relative mass difference between an 
initial target nucleus and a particular product nucleus for Ep — 1500 MeV. Letters are 
explained in Table 6. 
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